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ABSTRACT 

Context. Historically, due to instrumental limitations and a lack of disk detections, the structure of the transition from the envelope to 
the rotationally supported disk has been poorly studied. This is now possible with ALMA through observations of CO isotopologues 
and tracers of freezeout. Class 0 sources are ideal for such studies given their almost intact envelope and young disk. 

Aims. The structure of the disk-envelope interface of the prototypical Class 0 source, VLA1623A, which has a confirmed Keplerian 
disk, is constrained through modeling and analysis of ALMA observations of DCO^ (3-2) and C*®0 (2-1) rotational lines. 

Methods. The physical shucture of VLA1623 is obtained from the large-scale SED and dust continuum radiative transfer modeling. 
An analytic model using a simple network coupled with radial density and temperature profiles is used as input for a 2D line radiative 
transfer calculation for comparison with the ALMA Cycle 0 12-m array and Cycle 2 AC A observations of VL A1623. 

Results. The DCO^ emission shows a clumpy structure bordering VLA1623A’s Keplerian disk. This suggests a cold ring-like struc¬ 
ture at the disk-envelope interface. The radial position of the observed DCO^ peak is reproduced in our model only if the region’s 
temperature is between 11 K and 16 K, lower than expected from models constrained by continuum data and source SED. Altering 
the density profile has little effect on the DCO^ peak position, but increased density is needed to reproduce the observed C**0 tracing 
the disk. 

Conclusions. The observed DCO^ (3-2) emission around VLA1623A is the product of shadowing of the envelope by the disk observed 
in C**0. Disk-shadowing causes a drop in the gas temperature outside of the disk on >200 AU scales, encouraging the production 
of deuterated molecules. This indicates that the physical structure of the disk-envelope interface differs from the rest of the envelope, 
highlighting the drastic impact that the disk has on the envelope and temperature structure. The results presented here show that DCO^ 
is an excellent cold temperature tracer. 

Key words, stars: formation - stars: low-mass - stars: protostars - ISM: individual objects: VLA1623 - methods: observational - 
techniques: interferometric 


1. Introduction 


Rotationally supported disks have been observed extensively 
among most protostellar and pre-main sequence evolutionary 
stages ( Li et al.|2014| l. Recent studies have revealed the exis¬ 
tence of such disks in the Class 0 deeply embedded phase (Choi 


et al.|2010l|Tobin et al.|2013l|Murillo et al.|2013l|^della et al. 


2014 Lee et al.|2014 l. In the early stages of star formation, the 

envelope is not yet dispersed and contains enough material to 
influence the evolution of the star-disk system. The boundary 
between the disk and the envelope, known as the disk-envelope 
interface, must then play a role in the formation process. This 
region, however, is largely unexplored owing to limitations in 
the resolution and sensitivity, as well as the lack of observed ro¬ 
tationally supported disks, until now. Class 0 sources with con¬ 
firmed rotationally supported disks grant us the opportunity to 
study the chemical and physical structure of the disk-envelope 


* Current address: Universitat Heidelberg, Zentrum fiir Astronomie, 
Institut fur Theoretische Astrophysik, Albert-Ueberle-Str. 2,69120 Hei¬ 
delberg, Germany 


interface region, which is crucial for the next step of understand¬ 
ing this region’s role in star formation ( Sakai et al.|2014) l. 

Whilst CO isotopologues are good tracers of rotationally 
supported disks, they freeze out onto dust grains below the evap¬ 
oration temperature usually between 30-20 K pprgensen] 
et al.| 2005 i. These low temperatures are reached at the edge 


of the embedded disk ( |Visser et al.||2009) l. As a result molec 
ular species whose abundance is enhanced at low temperatures 
are needed to trace the disk-envelope interface. DCO^ emission 
is known to be optically thin and its abundance is enhanced at 
a narrow range of temperatures below the CO freeze-out tem- 
perature tracing the so-called CO snowline (e.g., Wootten|1987 


Roberts et al. 2003 


2013 I. Thus, DCO^ is a gooc 


Mathews et al. 

candidate molecule to trace the chemical, physical and kinematic 
structure of the disk-envelope interface. 

VLA1623-2417 (hereafter VLA1623) is a triple non-coeval 
protostellar system located in p Ophiuchus at ~ 120 pc ( |Murillo| 
& Lai|20T^. VLA1623A is the prototypical Class 0 source and 


emits predominantly in the (sub)millimeter range (Andre et al. 


|1993[ |Murillo & Lai|20r3| ). Lrom modeling of ALMA Cycle 0 
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Fig. 1. DCO^ compared with continuum and C’^O toward VLA1623A & B, marked with a square and triangle, respectively. The black ellipse 
shows the synthesized beam. Intensity integrated DCO^ (color-scale) with (a) 1.3mm continuum and (b) C**0. DCO^ peaks at 0.0433 Jy beam^' 
km s“’. Contours are in steps of 3, 5, 10, 15, 20, 40, 60 and 78cr with cr = 1 mjy bearn^* for 1.3 mm continuum, and -5, -3, 3, 5, 10, 15, 20 and 
25(T with (T = 13 mJy beam^' km s“' for C’*0. Velocity map (moment 1, color-scale) for (c) DCO^ and (d) C***0 overlaid with the corresponding 
intensity integrated map (contours). Contours are in steps of -10, -7, -4, 4, 7, 10 and 1 Icr with cr = 3 mJy beam^' km s“' for DCO^, and the same 
as in (b) for C'®0. 


C'^O observations. 


(20131 found that VLA1623A 


, Murillo et al. _ 

supports a Keplerian disk with a radius of at least 150 AU and 
a central mass M, of 0.2 Mq. The ALMA Cycle 0 observations 
also detected DCO^ molecular line emission toward VLA1623A 
bordering the C'^O disk. This grants us the opportunity to probe 
the disk-envelope interface of VLA1623A. 


In this paper we present the results of our ALMA observa¬ 
tions and simple chemical modeling, aiming to understand the 
physical structure of the boundary between the envelope and the 
disk in a Class 0 protostar. 


2. Observations 


Using the Atacama Large Millimeter/submillimeter Ar¬ 
ray (ALMA) we observed VLA1623 (pointing coordinates 
a=16;26:26.419 5=-24:24:29.988 J2000) during the early 
science Cycle 0 period on 8 April 2012. Observations were 
carried out using the extended configuration, comprised of 16 
antennas with a maximum baseline of ~400 m, in Band 6 (230 
GHz). Total observing time was one hour with a 58% duty cycle. 
The spectral configuration was set to observe DCO^ (3-2), C^^O 
(2-1), and '^CO (2-1) in addition to continuum. N 2 D^ was also 
observed but no significant detection was made. Data calibration 
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Fig. 2. Intensity (contours) and velocity (color-scale) integrated maps 
of DCO^ ACA observations. Note the 10 times larger scale of this fig¬ 
ure compared with Fig.[^ The positions of VLA1623A, B and W are 
marked with a square, triangle and circle, respectively. The size of the 
filled square shows the size of the C'*0 disk. The red and blue lines in¬ 
dicate the outflow direction while the gray line shows the mayor axis of 
the disk plane. The black ellipse shows the synthesized beam. Contours 
are in steps of -3, 3, 5, 7, 10, 15 and 17cr with cr = 220 mjy bearn^* 
km s“'. 
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velocity (km/s) 

Fig. 3. Cycle 0 DCO^ (black) and C**0 (gray) PV diagram. Cuts are 
made at P.A. = 26° for DCO^ and P.A. = 35° for C‘*0. Keplerian (solid 
line) and Infall plus Keplerian (dashed line) rotation PV models with 
M, = 0.2 Mq are overplotted l |Murillo et al.|2013) . Contours are in steps 
of -3, 3, 5, 10 and 15cr with a = 19 mJy beam“‘ for C**0, -7, -5, -3, -2, 
2, 3, 4, 5cr with cr = 12 mjy bearn^* for Cycle 0 DCO^ 


was done with 3C 279, 1733-130, and Titan for bandpass, gain, 
and flux calibration, respectively. 

In this paper, we present the results and analysis of the DCO^ 
(3-2) (rest frequency: 216.11258 GHz) observations. The DCO^ 
data were calibrated jointly with the continuum, C'^O (2-1), and 
*^CO (2-1) data. Further calibration details and results from the 
other observed lines can be found in |Murillo et al.| ( [2013| ). The 
spectral set-up provided a velocity resolution of 0.0847 km s '. 
The synthesized beam size for the DCO'^ images is 0.85" x 
0.65" with P.A. = 96°. The rms noise of the channel map is 12 
mJy beam'* for a spectral resolution of 0.0847 km s'*, giving a 
peak S/N = 7. 

Our ALMA observations provide a maximum scale of 4" and 
a held of view (FOV) of 24", with emission between 4" and 24" 
largely Altered out by the interferometer. The FOV together with 
the beam size of 0.85" constrain the scale to which any analysis 
of the data can be done. 

In addition to the ALMA Cycle 0 observations, we present 
the DCO^ (3-2) results from our ALMA Cycle 2 Atacama Com¬ 
pact Array (ACA) observations carried out on 7 August 2014 
(pointing coordinates a=16:26:26.390 d=-24:24:30.688). Total 
observing time was 2 hours. Data calibration was done with 
J1517-243 and Mars for flux, 11625-2527 and Mars for gain, and 
11733-1304 for bandpass. The rms noise of the DCO^ channel 
map is 73 mly beam'* for a spectral resolution of 0.021 km s'*, 
with a synthesized beam of 8.6" x 4.2" with P.A. = -76°. These 
observations, with a mosaicked area of 6', provide the DCO*^ 
emission between 4" and 18" scales. 


3. Results 

3.1. ALMA 12-m Array Cycle 0 

The detected DCO^ emission (Fig. is located between a ve¬ 
locity range of 2.8 to 5.2 km s'* and shows a clumpy structure 
with two main clumps to the north and south of VLA1623A. 
The southern clump emission is stronger and is offset by about 
2.5" (300 AU) from VLA1623A and in addition it borders the 
red-shifted emission of the disk traced in C*^0. The northern 
clump slightly overlaps the blue-shifted emission of the C*^0 
disk and borders VLA1623B’s continuum emission. A couple of 
clumps with emission between 3 and lOcr are observed near the 
continuum peaks of VLA1623A & B. No significant emission 
was detected toward VLA1623W, separated by 10" to the west 
from VLA1623A, possibly either because of the lack of DCO*^ 
or because the emission is too weak and Altered out. 

The C*^0 line emission tracing the disk was found to be in¬ 
fluenced by the outer envelope, with the blue-shifted emission 
being affected more than the red-shifted emission ( Murillo et al.| 
2013|l. We expect the same to hold true for the DCO^ emission, 
thus potentially explaining why the northern blue-shifted clump 
emission is weaker than that for the southern red-shifted clump 
(Fig.Qa&b). 

The velocity weighted (moment 1) map of the DCO^ emis¬ 
sion (Fig. [^) shows that its velocity gradient is similar to that of 
C*^0, with the northern clump being blue-shifted and the south¬ 
ern clump being red-shifted, but with a smaller velocity range. 
The Position-Velocity (PV) diagrams of DCO*^ and C**0 emis¬ 
sion (Fig.[^ are constructed and over-plotted with the best fitting 
thin disk models, Keplerian and Infall plus Keplerian out to 150 
AU, obtained by |Murillo et al.| ( |2013| l. It appears as though both 
line emissions are well described by pure Keplerian rotation out 
to 300 AU. However the DCO^ emission is too weak, peaking at 
Icr in the channel map, to carry out further kinematical analysis. 
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Fig. 4. Cartoon of the chemical network used to analyze the observed 
DCO^ emission. Red denotes the target molecule. Green the input 
molecular concentrations and orange the species involved in the bot¬ 
tleneck reaction. 


The similar velocity signatures of the emission from both species 
may be due to the disk edge dragging along material from the en¬ 
velope. 

The above simple analysis of the observations indicates that 
the detected DCO^ emission may be tracing a ring which borders 
the C'*0 disk. In addition, DCOM’s velocity gradient and PV di¬ 
agram suggest that the ring is undergoing Keplerian rotation, this 
is even more plausible given that the (2-1) emission traces 
a rotationally supported disk with a Keplerian velocity profile 
HMurillo et al.|20131 l. 


3.2. ALMA ACA Cycle 2 

The DCO^ (3-2) emission traced with the ACA confirms that the 
emission is concentrated around VLA1623 (Fig. |^, it does not 
peak on VLA1623A and instead encircles it, forming a shell-like 
structure. The ACA detects only weak emission at the position 
ofVLA1623W. 

The DCO^ emission mapped with the ACA shows the same 
velocity gradient (Fig. as on the small scales, indicating that 
the kinematic structure is the same throughout. 


4. Analysis 

DCO^ is detected at the disk-envelope interface of VLA1623A. 
Given that DCO^ is optically thin and a good probe of tempera¬ 
ture and CO freeze-out regions, we model the observed emission 
aiming to probe the physical and chemical structure of the disk- 
envelope interface. In this section we describe the model used 
and the results obtained from the model. 
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4. 1. DCO* chemical network and model 

We model the observed DCO^ emission using a simple chemical 
network. Such a network, while it may not account for every pos¬ 
sible reaction, is preferred as it gives insight into how the physi¬ 
cal parameters of the temperature and density affect the observed 
emission. The network used is a steady-state analytic model that 
only takes the basic reactions that lead to DCO^ production and 
destruction into account (Table[^ Fig.|^. Since DCO^ is formed 
by the reaction of and CO, the rate-determining reactions 
in our network are 

H+-HHD—>H2D+-hH2 (1) 

with the back reaction 


H 2 D+ -H H 2 -H AE —> 4- HD (2) 

where the activation energy due to the difference in zero point 
energy is AE ~ 220 K. Eor reaction 5 in our network (Table [^, 
for simplicity, we adopt the total rate coefficient summing over 
all three pathways of the reaction. The rate coefficients for a two- 
body reaction are given by 

where T is the gas temperature, while 

k = ^s-' (4) 


gives the rate coefficient for cosmic-ray ionization (reaction 1 
from Table [T]). 

Our model takes as input a source density and temperature 
profile as a function of radius and the parameters needed to cal¬ 
culate the rate coefficients. The CO evaporation temperature T^y, 
desorption density nje and CO abundance Xqq are free parame¬ 
ters. The CO abundance is assumed and not calculated. The CO 
evaporation temperature dictates when CO is in the gas phase 
(T > Tev) or freezes onto the dust grains {T < T^y). In a simi¬ 
lar manner, the desorption density sets the boundary when the 
freeze-out timescales are too long (n < Ude) compared to the life¬ 
time of the core pprgensen et al.|2005| l. We assume the density 
profile is equal to the H 2 density mh, and the abundance of HD 
2fHD = 10^^ with respect to the total hydrogen nuclei density rifj 
- 2n(li2). The model returns the calculated concentrations as a 
function of radius. The results can then be input into excitation 
radiative transfer programs such as RATRAN ( [Hogerheijde &| 
van der Tak|2000 1 for further analysis. 

Given that CO is one of the parent molecules of DCO^, we 
study the effect of CO abundance through the use of different 
abundance profiles following the models detailed in J0rgensen| 
et al. (20051 and |Yildiz et al.| ( 2010| l. The possible CO profiles 


are Constant and Drop abundance (Fig.|^. The former represents 
a fixed CO abundance throughout the core, whereas the latter, 
constrained by T^y and nje, is used to account for CO freeze-out 
in the chemical network. Eor the constant profile, the abundance 
is denoted by Xq. Eor the drop profile, the abundance in the inner, 
drop and outer regions are denoted by Xi„, Xu and Xq, respec¬ 
tively. Previous studies using multiline single dish C'^O have 
found that the abundance at Ain is lower than Aq for a number of 
sources ( [Alonso- Albi et al.|2010[rYildiz et al.|2010||2013| l. One 
explanation is that some fraction of the CO ice is transformed 
into more complex and less volatile carbonaceous species in the 
cold phase. We thus take this effect in our model into account. A 
lower abundance of CO, due to freeze-out, allows an increase in 
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Table 1. DCO^ chemical network reactions and adopted rate coefficients 


ID Reaction 

r 



/ 

Ref. 

1 n ^ :: . fT+ . 

s-' 

3 -1 

cm s 


K 



2 

3 

4 
4a 
4b 

5 

6 

7 

8 

9 


H+ -H H, - 
-hHD 


m 


+ H 




-hH, 




-hH, 


H, 


H2D^ -H P-1X2 
+ 0 -H 2 


H 


-I- e - 
3 -H CO — 
HCO+ H- e - 
H 2 D+ -H CO 
DCO+ -I- e - 


-i-HD 

^ H+ +HD 
^ H+ + HD 

^ HD-hH 
HCO+ -H H2 
H-hCO 
DCO+ H- H2 
D-i-CO 


2.08 (-9) 
3.50 (-10) 
3.50 (-9) 

1.40 (-10) 
7.00 (-11) 

6.79 (- 8 ) 
1.61 (-9) 

2.80 (-7) 
5.37 (-10) 

2.40 (-7) 


0 

0 

0 

0 

0 

-0.52 

0 

-0.69 

0 

-0.69 


0 

0 

220.0 

232.0 

61.5 

0 

0 

0 

0 

0 


2 

3,4 

3,4 

5 

5 

5" 

6,7 

8 

4 

4 


Notes. The reactions in bold (4a and 4b) substitute the back reaction of the bottleneck (4) in the network when o-Hj and p-Hj are included. 

^ is the cosmic-ray ionization rate of Hj. a is the rate coefficient at T = 300 K. jj is the temperature coefficient. 7 is the activation 
barrier. Value obtained by summing the rate coefficients for the reaction channels to three possible products. 

References. (l)|Black|1975| (2)|Theard & Huntress| 19741 (3)|Gerlich et al.|20021 (4)|Albertsson et al.|20131 (5)|Walmsley et al.|20041 ( 6 )|Plessis| 
let al.|2010| (7) OSU Database; (8) |Amano|19'^ 


the abundance of H 2 D^ (Mathews et al.|20T3i. However, since 


both molecules, CO and H 2 D^, are parent molecules of DCO 
balance must be reached before the effective formation of DCO^ 
takes place. This scenario is found to be common for the en¬ 
velopes of early embedded protostars pprgensen et al.|[2005| l 
where the outer region shielded from the protostellar and in¬ 
terstellar radiation heating has low enough temperature that CO 
freezes out onto dust grains. 

In molecular clouds, and consequently in protostellar cores, 
H 2 chemistry plays a major role, hence the oitho-to-para ratio of 
H 2 influences the chemical reactions, and has been found to be 
crucial to the deuterium chemistry (|Flower et al.||2006} |Pagani| 


et al.||20()9 1 . The effect of ortho- and para-H 2 ( 0 -H 2 and P-H 2 ) y, 


IS studied in our model through the inclusion of the ortho-to- 
para ratio (o/p) and the distinction of 0 -H 2 and P-H 2 in the back 
reaction of the bottleneck (Eq. ^ in the chemical network, o- 
H 2 and P-H 2 are only added in the back reaction since it is here 
where the distinction has a significant effect (Table [^l. We set a 
lower limit on o/p of 10 “^ at low temperatures, as constrained by 
models and observations ( jFlower et al.|200^|Faure et al.|2013l l. 
The 0 -H 2 and P-H 2 reactions and their parameters for the rate 


coefficients are taken from |Walmsley et ak ( 2004| l. 

When 0 -H 2 and P-H 2 are included in the network, a thermal 
(LTE), upper- or lower-limit o/p ratio can be selected. In LTE, 
the ortho-to-para ratio is given by 


o/p = 9 exp 


-170 


( 5 ) 


where T is the gas temperature. Selecting the upper-limit ratio 
produces 3 times more 0 -H 2 than P-H 2 . Since the back reaction 
with o-Hj has a lower activation barrier 7 than with P-H 2 , the o/p 
upper-limit implies that H 2 D^ is being destroyed faster tha n gen - 
erated, leading to a decreased production of DCO^ (Eig. A.3 1 , 


since H 2 D" is a parent molecule of DCO^ (Eig. ffl. The lower- 
limit ratio, on the other hand, implies more p-Hj which has a 
higher activation barrier 7 for the back reaction, thus H 2 D^ is 
generated faster than it is destroyed in turn increasing the DCO^ 
production (Eig. |A.3| . 

As a starting point for our analysis of VLA1623A’s DCO^ 
emission, we use the density and temperature profile of 
VLA1623 obtained by |j 0 rgensen et al^ ( |2002| l where 30 K is at 



Radius (AU) 

Fig. 5. CO abundance profiles used in the model. The vertical dashed 
lines show the limits for the Drop abundance profile, evaporation tem¬ 
perature Tev and desorption density Xq denotes the abundance in 
the constant profile. Vin, Xp and Aq are the inner, drop and outer region 
abundances for the drop profile. 


~1.5" assuming a distance of 120 pc. The profile was obtained by 
fitting single dish JCMT continuum images and the spectral en¬ 
ergy distribution (SED) with continuum radiative transfer mod¬ 
eling using DUSTY resulting in a power law density profile of 
the form n oc In the single dish continuum observations 

VLA1623A and B are unresolved and the density profile ex¬ 
tends well beyond the ALMA field of view. It is expected, how¬ 
ever, that VLA1623A dominates at 870 /tm and that VLA1623B 
does not contribute much to the 450 jum continuum ( |Murillo &] 
Lai||2013|l. Thus, th e density and temperature profile obtained 
by |j 0 rgensen et al^ ( |2002 1 is representative of VLA1623A since 
VLA1623B is not significantly contributing to the continuum 
emission or the SED used to constrain the profiles. 

The results of the analytic chemical network are run through 
the molecular excitation and line radiative transfer program RA- 
TRAN to generate line emission maps. Since the structure we 
are trying to reproduce is ring-like, we calculate the level popu- 
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Table 2. Input CO abundance profiles 


Case 

CO 

Lev 

«de 


Yd 

Yo 


profile 

(K) 

(cm“^) 




Fixed 

Constant 





10 -^ 

a 

Drop 

30 

X 

o 

10^5 

10-5 

10-4 

b 

Drop 

35 

3x10^ 

10-5 

10-5 

10-4 

c 

Drop 

35 

3x10'’ 

10-5 

5x10-5 

10-4 

d 

Drop 

35 

3x10'’ 

10-5 

10 -^ 

10-4 

e 

Drop 

35 

5x10^ 

10-5 

10-5 

10-4 

f 

Drop 

35 

3x10'’ 

10-5 

10 -*^ 

10-4 


Notes. See Fig.j^for definition of Xi,,, Xj^ and Xq 


lations with the 1 -D version of RATRAN, and then run the level 
populations with the 2-D ray tracing to form the ring structure. 
The produced spectral image cubes are convolved with a Gaus¬ 
sian beam with the dimensions of the synthesized beam, contin¬ 
uum subtracted and then an intensity integrated map is gener¬ 
ated. Radial profiles are extracted from the resulting images and 
compared with the observed profiles, which are integrated over 
the extent of the detected emission. 


We find no significant difference in the peak position be¬ 
tween the method used here and running the models through the 
ALMA simulator with the actual ALMA configuration or our 
Cycle 0 observations. However, the ALMA simulator does show 
that the emission >4" is indeed filtered out by the Cycle 0 obser¬ 
vations, producing somewhat narrower radial profiles. M olecular 
data for DCO^ and C'^O are obtained from LAMDA |Schoier 


et al. 


2005 


Flowep 

focuses 


5 on the 


C'^O: 


2010 


DCO extrapolated from 


Yang et al. 

. As DCO^ is optically thin, the comparison mainly 
tie position of the peak and the integrated intensity 


profile with respect to radius, the velocity profile has no effect 
on the resulting model integrated intensity. Thus, we assume a 
free-fall velocity profile with a central mass of 0.2 Mq ([Murillo | 


et al. 20131 and Ldust = Tgas. The RATRAN output maps are con- 
volved with the observed beam (0.85" x 0.65", RA. 96.24°) and 
compared with the observed emission through radial cuts. In ad¬ 
dition to comparing DCO^, we also compare C^^O in order to 
further constrain the physical structure of the disk-envelope in¬ 
terface. However, for C'^O we assume a pure Keplerian veloc¬ 
ity profile, in agreement with the rotationally supported disk it 
traces. For both emission lines, we use an inclination of 55° (90° 
= face-on), in accordance with the results obtained by |Murillo| 

|erar] ( |20T^ . 


4.2. Modeling results 

In this section the results of altering the CO abundance, density 
and temperature profiles in the chemical model are discussed. 
Figure [^presents the models, with the light blue region showing 
the location of the observed DCO^ peak emission with respect 
to VLA1623A’s position. Model naming follows the scheme xYz 
where x is the test number from Table Y is either C for con¬ 
stant or D for drop abundance, and z is the case from Table 


4.2.1. Chemical properties 

To probe the chemical conditions, we establish one case for the 
constant CO abundance with Xqo - 10 and six cases for the 
drop CO abundance with Ad ranging from 5x10 ® to 10“* and 
varying T„ and Wde- Parameter ranges were selected based on 
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trends found in previous work pprgensen et al.|2005[|Yildiz et al.| 
2010|l and adapted to the current observations. Parameters for 


each case are listed in Table As a zeroth order test, we com¬ 
pare the abundance profiles to the radial position of the observed 
DCO^ peak in the ALMA 12-m array data in order to investi¬ 
gate which chemical conditions best approximate our observa¬ 
tions. This comparison neglects the fact that the peak emission 
radius also depends on the DCO^ excitation, which will be taken 
in Section l43] into account. 

Figure [^presents the model abundance profiles as functions 
of radius for different assumed CO abundances and physical 
structures. In general, the concentrations of CO and HCO^ drop 
with radius whereas those of H 2 D^ and DCO^ increase with ra¬ 
dius due to the lower temperature farther away from the source. 
Constant CO abundance (Fig. left column) does not appear to 
produce a DCO^ peak within the expected region. Since the ra¬ 
dial abundance profile of DCO^ is not altered by changing Aco, 
we focus instead on the drop CO abundance profile. 

While the drop CO abundance profile cannot alone alter the 
position of the peak, it produces several trends interesting to 
note. A kink in the abundance of DCO^ forms at for all the 
cases examined, though in most cases it is relatively small. Al¬ 
tering the abundance in the drop Ad changes the shape of the 
peak but d oes not significantly alter its position nor its abundance 
(Fig. A.2 1 . Decreasing the Ad below 10“^, however, causes the 
abundance of the peak to drop by several orders of magnitude 
and become similar in magnitude to the kink at producing 
two peaks which are not observed (Fig, |A.2[ case f). Varying Tev 
and rids, in order to constrain the width of the drop, generates a 
very narrow peak which increases and drops quickly. 

To test the effect of the o/p ratio on our chemical network, we 
vary the ratio from the thermalized value to the upper- and lower- 
limit for the drop CO abundance profile. In all cases, we find that 
setting o/p to the lower-limit does not change the position of the 
DCO^ emission since the peak and bulk of the DCO^ concen¬ 
tration for t he th ermalized o/p is already located in the o/p=10“^ 
range (Fig. A.3 i. In fact, the peak of the modeled DCO^ emis¬ 
sion starts to decrease as the o/p ratio increases. The lower-limit 
o/p only alters the inner regions (<100 AU), which the data do 
not constrain. The upper-limit reduces the overall production of 
DCO^ and pushes the peak outward to larger radii. In conclu¬ 
sion, altering the CO abundance and o/p ratio does not produce 
the observed results. 


4.2.2. Physical properties 

Altering the chemistry of the model does not reproduce the ob¬ 
served DCO^ peak position, regardless of the case used. Thus we 
alter the density and temperature profile in order to find the con¬ 
ditions necessary to reproduce the observed emission. We set up 
4 tests, including the original profile from |j0rgensen et ar] ( |2002] l 
referred to as Test 1, by increasing or decreasing by a constant 
factor either the density or temperature profile. Parameters for 
each test are listed in Table |3] 

Increasing or decreasing the density profile by one order of 
magnitude. Tests 2 and 3, does not generate a significant ra¬ 
dial shift of the DCO^ peak compared with the original pro¬ 
file (Fig. 1^ middle row). Altering the density, however, affects 
the concentration of the peak upwards or downwards with an in¬ 
crease or decrease of the density, respectively. 

The factor of 1.5 decrease in the temperature profile with 
cases b and d with Ad 10 ® - 10“^ shifts the modeled DCO^ 
peak effortlessly to the observed position (Fig.[^. Decreasing the 
temperature profile between a factor of 1 to 2 moves the modeled 
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Fig. 6. Chemical network model results. Left and center columns show the results for constant and drop CO abundance, respectively. Top, middle 
and bottom row show the results of Test 1, 2 and 4, respectively. Light blue shaded region shows the observed peak position of DCO^. Right 
column shows the input source profile, where the original DUSTY source profile is shown with a dashed line so as to compare with the input for 
the specific test. 

Table 3. Input source profiles for each test 


Test 

7’4au (K) 

n4Au (cm^-^) 

Source Profile 

1 

250.0 

1.62 X 10'^ 

from J0rgensen et al. (|2002|) 

2 

250.0 

1.62 X 10'° 

Density increased 1 order of magnitude 

3 

250.0 

1.62 X lO*' 

Density decreased 1 order of magnitude 

4 

166.7 

1.62 X 10° 

Temperature decreased by a factor of 1.5 with Tiowiim = 8 K 
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Fig. 7. Comparison between observations and radiative transfer model¬ 
ing of the chemical network results. Here the best approximations are 
highlighted, with model 4Dd being the preferred model. The models are 
referred to by test, abundance profile (Constant “C" or Drop “D") and 
case. Thus, 4Dd is the model for test 4 with Drop abundance and case 
d. Black line shows the DCO^ observations integrated over the southern 
clump. Gray solid and dashed lines show the Icr and 4cr levels, respec¬ 
tively. 


DCO peak inward (Fig A.4 1 . For the constant CO abundance, 
any alteration of the temperature either over-estimates the mod¬ 
eled DCO^ abundance or produces a peak too far inward. For the 
drop CO abundance, a temperature drop less than a factor of 1.5 
does not move the peak inward enough, whereas a larger factor 
moves it too far inward. Cases b and d with Xo =» 10“® - 10^^ 
shifts the modeled DCO^ peak effortlessly to the observed po¬ 
sition (Fig. 6 bottom row). Changing the CO abundance in the 
drop (Table 2 1 does not alter the location of the peak, as expected 
from the results presented in Sec. 4.2. 1| Examining the temper¬ 
ature profile shows that the DCO^ emission peaks at a range of 
11-16 K. This is lower than the expected 20 K at ~3" inferred 
from radiative transfer modeling of the observed continuum data 
Pprgensen et al.|[2002| l. We limit the decrease in the tempera¬ 
ture profile to not fall below 8 K ( Zucconi et al.|200T I, however 
this limit generates no significant change in the outcome of our 
model since the limit falls near and beyond the edge of the FOV 
of our observations (Fig.[^ bottom right panel). 


From the results obtained from altering the chemical and 
physical conditions in our model, we can deduce that the DCO^ 
emission is located at its observed position due to physical condi¬ 
tions, namely a lowered temperature, and not to special chemical 
conditions in the disk-envelope interface. 


4.2.3. Comparison with full chemical network 

To explore the limitations of the simple network, we compare 
the results to a full time-dependent deuterated chemical net¬ 
work, based on the RATE06 version of the UMIST Database for 
Astrochemistry ( [Woodall et al.||2007|l extend ed with deuterium 
fractionation reactions ( [McElroy et al.|2013| l. Models with only 
gas phase chemistry and including gas-grain balance (freeze- 
out, thermal desorption, and cosmic-ray-induced photodesorp¬ 
tion) are run. The models start with the same initial abundan- 
dances (HjiHDiCO = I;3xl0“^:lxl0“^) and the same physical 
structure as the simple network. We find that the large network 
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Fig. 8 . Same as in Fig. ^ but for C**0, with 2Dd being the preferred 
model. Gray solid and dashed lines show the Icr and 3cr levels, respec¬ 
tively. 2Dd, increased density profile, provides the best approximation 
to the observed C**0. 


confirms both the general trend found in the simple network and 
also reproduces the peak position and abundance of DCO^ at the 
observed position for abundances extracted at early times (~ 10 ^ 
yr). 


Inspection of the full network shows that in addition to the 
reactions listed in Table the HCO^ + D —> H + DCO^ reac¬ 
tion can become important if D/HD is large. The back reaction 
has a reac tion barrier of ~800 K, thus leading to fractionation 
of DCO^ ( Adams & Smith|[l985 i. In practice, this reaction is 
only relevant for much lower densities than encountered in the 
VLA1623A envelope model. Similarly, dissociative recombina¬ 
tion of Hj with electrons cannot be neglected at low densities. 
At very early times (< 10^ yr, depending on density), reactions 
with and D^ become more significant than those with H 3 
and H 2 D^. Finally, in full gas-grain models the bulk of the CO 
is frozen out in the cold outer envelope resulting in very low 
DCO^ concentrations, unless an efficient non-thermal CO des¬ 
orption process is included. 


4.3. Emission prof lies 

Aiming to better constrain the physical conditions of the disk- 
envelope interface we compare the models and observations of 
both DCO^ and C'^O emission computed with RATRAN. For 
the comparison, the observed southern red-shifted clump in the 
ALMA 12-m array data for both lines is selected. The reason is 
two-fold; i) the red-shifted emission of the C'^O disk suffers less 
absorption from the outer envelope than the blue-shifted lobe as 
noted above ( [Murillo et al.|[20T3] l; ii) the southern red-shifted 
clump traced by DCO"^ is the strongest and most prominent. 
Similarly, for the model of C'^O only one lobe of the disk is 
selected for comparison. 

Consistent with the abundance plots, the constant CO abun¬ 
dance profiles do not fit the observed DCO^ emission in any of 
the four tests, producing very broad peaks between 3" and 5" 
away from the source for tests 1, 2 and 3; and a peak at 1" for 
test 4. In a similar manner, for the drop CO abundance scenario, 
tests 1, 2 and 3 all produce peaks beyond 3" with varying broad¬ 
ness, and thus do not approximate the observed DCO^ emission. 
This occurs for all examined cases (Fig. ig. Models 4Db and 
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4Dd approximate the observed DCO^ emission well (Fig. [t] or¬ 
ange and red lines, respectively). These results are in agreement 
with the conclusions drawn from the concentrations in the ana¬ 
lytics chemical network model (See Sec. |4.2.2| ). 

Comparing the results of the model with the observed C'^O 
emission (Fig. [^, we find that in all four tests a constant CO 
abundance over-predicts the extent of C**0, with the emission 
well above our her level extending out to 3" or further, whereas 
we observe it only out to less than 2" from the central source. 
When the drop CO abundance is introdu ced, the source profile of 


VLA1623 from Jprgensen et al. (20021 overestimates the C O 
emission in all cases, and generates a second peak at about 3" 
at half the intensity of the central peak in almost all cases. The 
width of the drop does not have any significant effect on the mod¬ 
eled emission. The decreased density and decreased temperature 
profile tests do not fare better, as they again largely overestimate 
the C'^O beyond 2" and even produce secondary peaks. Case a 
(Ad = 10 “^) for the original profile and decreased density seem 
to show some promise, however the modeled emission is around 
3cr at 1.5" whereas observations at that radius are closer to Icr. 
Surprisingly, we find that the increased density profile produces 
the best results with case d (Aq = 10 “^), the same as that for 
DCO+ (Fig.§ red line). 

Hence we find that the DCO^ and C'^O observed emission 
are not reproduced by the same physical structure (Figs. [7] and 
[^. DCO^ is well modeled by a temperature profile a factor of 
1.5 lower than that needed to model C'^O, whereas C^^O is re¬ 
produced by a density profile an order of magnitude higher than 
required for DCO^. The abundance in the drop Ad necessary to 
reproduce the emission is the same for both lines. Ad=10 ^^. 


4.4. Large scale emission vs. model 

We also compare the radial integrated intensity profile of the 
DCO^ ACA observations with the model, both along the plane 
of the disk and perpendicular to the disk in the outflow direction, 
since these data probe larger spatial scales than the 12-m array 
data. For this purpose, the models were run through RATRAN’s 
ray tracer sky again but instead using the cell size of the ACA 
observations (0.84") and then were convolved to match the reso¬ 
lution of the ACA observations. The observed radial profiles are 
extracted from the red-shifted emission. 

The profile along the plane of the disk (peak at ~4") is well 
described by the same model as in the small scale, the drop CO 
abundance profile for case d, with the temperature decreased by 
a factor of 1.5 (4Dd, Fig. The temperature decreased by a 
factor <1.5 does not provide a better estimate. The constant CO 
abundance again overestimates the amount of DCO^ produced 
in the outer regions. 

On the other hand, the integrated intensity profile along the 
outflow direction (peak at distances =s5") is very well described 
by the model 1C (Fig.[^. This indicates that in regions that are 
not shadowed by the disk, the production of DCO^ is as expected 
from dust continuum models. 


5. Discussion 

The results of our modeling show that the position of the DCO^ 
emission (~2.5" = ~300 AU) is closer to VLA1623A than ex¬ 
pected based on the spherically symmetric dust continuum ra¬ 
diative transfer modeling from DUSTY (~4" = ~480 AU). The 
results further show that the observed emission peaks at a dust 
temperature range of 11-16 K. Comparison of the radiative trans¬ 




Fig. 9. Radial profiles of the DCO^ ACA observations along the disk 
plane (top) and along the outflow axis perpendicular to the disk (bot¬ 
tom) as shown in Fig. Radial profile of the DCO^ Cycle 0 observa¬ 
tions is shown in light gray. Observed radial profiles are overlaid with 
models convolved to the resolution of the ACA observations. Gray solid 
and dashed lines show the Itr and 3cr levels, respectively, of the ACA 
emission. These results evidence that the DCO^ emission along the disk 
plane is best approximated by model 4Dd, whereas along the outflow is 
best described by model 1C. 


fer results from RATRAN to the C^*0 emission from the disk 
suggest that the disk has a density higher by one order of magni¬ 
tude than the emitting structure traced by DCO^. 


A possible explanation would be that VLA1623A is a Very 
Low Luminosity Object (VeLLO, [Young et al.||2004l |Dunham| 
jet al.|2008] l undergoing episodic accretion and just coming out 
of the quiescent phase. This is highly unlikely, however, since 
[Johnstone et al.j P013[ l find that the timescale for dust and gas 
to heat up after an accretion burst is short, on the order of hours 
to weeks. VLA1623A has been reported of having a bolometric 
luminosity between 0.4-2 Lq from early observations to more 
recent work ([Andre et al.|199^ |Froebrich|2005[ [Murillo & Lai| 


20T3)[Chen et al.[[2013[ ), in contrast to the expected luminos¬ 

ity of 10“^ Lq for VeLLOs. Hence VLA1623A has not recently 
come out of the quiescent accretion phase and the location of the 
observed DCO^ emission can not be attributed to the relic of a 
previous phase of decreased accretion. 
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Fig. 10. Radiative 2D models based on |Harsono et al. i2013 i. For all 
models an envelope with M^m = 1 Mq irradiated by a 1 Lq protostar is 
assumed, having an outflow cavity with a full opening angle of 30°. The 
dashed lines indicate the temperature contours for 100, 50 and 20 K. The 
top and middle panels include a 180 AU and 0.02 Mq disk assuming a 
0.1 AU scale height, with centrifugal radius = 200 AU (top) and 50 
AU (center) orientated in the horizontal direction. The outflow cavity is 
oriented in the vertical direction. The bottom panel shows the envelope 
without a disk. 


A more plausible explanation to the position and tempera¬ 
ture of the region containing the DCO^ emission is disk shield¬ 
ing. DCO^ is observed to border VLA1623A’s C'^O disk in our 
present data. This could shield the outer parts of the disk from 
heating by the central protostar, thus moving in the freeze-out 
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zone of CO, enhancing the production of DCO^ closer to the 
protostar. This scenario is further supported by the result of our 


simple chemical model that the C'^O disk is more dense than the 


region of the DCO emission and than expected from Jprgensen 
et al. (|2002|l’s envelope density profile of VLA1623. 


We test the possibility of the disk-shielding scenario with 
radiative 2D disk plus envelope models using radiative transfer 
methods as in Harsono et al. (2013[l (Fig. 101. For the models, we 


assume a central protostellar source of 1 Lq (Murillo & Lai 2013 
Chen et al.|2013 i, an envelope with a mass of 1 Mq (Andr6 et al.| 


1993[ [Froebrich 2005| l, with the addition of a disk mass of 0.02 
Mq and radius of 180 AU ( [Murillo et al.|2013] l. The outflow cav¬ 
ity is assumed to have an opening angle of 30°. A thin disk with 
a scale height of 0.1 AU is adopted. Disk flaring is not included 
since we have no information on the flaring of VLA1623A’s disk 
and the thin disk model approximates the C**0 kinematics well. 
Two values for the centrifugal radius 200 AU and 50 AU, are 
chosen. Within the centrifugal radius the velocity structure of 
the disk is Keplerian. We find that even for such a thin disk the 
temperature along the plane of the disk is lower than for the en¬ 
velope at the same radius (Fig. top and middle panels), thus 
moving the CO freeze-out zone closer to the protostar along the 
edge of the disk than in other regions of the core. For either cen¬ 
trifugal radius the temperature beyond 200 AU drops well below 
20 K. Finally, we test whether the presence of a disk makes a 
difference in the location of the CO freeze out region. Figure [T0| 
bottom panel, shows the model with the same conditions as in 
Figure [TO middle panel, but without the disk. This shows that 
omitting the disk causes the CO freeze out region to move out¬ 
ward to ~400 AU, 150 AU further out than the models including 
a disk. 


The temperature conditions required to reproduce the ob¬ 
served DCO^ emission with our simple chemical model are 
therefore in agreement with the results obtained from the 2D ra¬ 
diative transfer disk plus envelope model. These results strongly 
support the scenario where a disk can shield the regions at its 
edge from heating by the protostar. This shielding causes the CO 
freeze-out region to move inward toward the edge of the disk, 
bringing low-temperature enhanced molecules, such as DCO^, 
closer to the central protostar along the plane of the disk. How¬ 
ever, the rest of the envelope is largely unaffected by the disk 
and a shell of molecules such as DCO^ forms at a radius pre¬ 
dicted by spherically symmetric models. The ACA observations 
of DCO^ toward VLA1623A provide further evidence for this 
effect. These results show that the DCO^ emission peaks closer 
to the source along the disk plane than along the outflow direc¬ 
tion. 


6. Conclusions 

This work presents the results and analysis of ALMA Cycle 
0 Early Science Band 6 observations of DCO^ (3-2) in the 
extended configuration toward VLA1623A probing subarcsec 
scales, as well as Cycle 2 ACA data probing the larger scales. 
A simple chemical network was setup taking ortho- and para- 
H 2 into account in the rate-determining reactions. The density 
and temperature profile of VLA1623A was obtained from fitting 
the SED and dust continuum data with radiative transfer mod¬ 
eling using DUSTY ( jjprgensen et al.||2002| l. Our simple chem¬ 
ical model coupled with VLA1623A’s physical structure served 
as input for line radiative transfer calculations with RATRAN 
with 2-D ray-tracing. The CO abundance, density and tempera¬ 
ture profiles were altered to study the effect of each parameter 
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Fig. 11. Cartoon of the disk and envelope of VLA1623A. The ring 
and disk, traced in DCO^ and C’*0 respectively, are observed with 
ALMA 12-m array, while the outer spherically symmetric shell, traced 
in DCO^, is observed with ALMA AC A. 


on the location of the observed DCO^ peak. The results of our 
observations and analysis can be thus summarized; 

1 . DCO^ is observed to border the C'^O disk around 
VLA1623A. Both emission lines show similar velocity gradi¬ 
ents (blue-shifted to the north and red-shifted to the south), with 
DCO^ emission between 2.8 to 5.2 km s '. The PV diagrams of 
C'^O and DCO^ suggest that both emission lines are well de¬ 
scribed by Keplerian rotation. However, the DCO^ emission is 
weak, thus no further kinematical analysis was carried out, and 
whether the disk extends out to 300 AU cannot be confirmed. 

2. Using a simple chemical network with the inclusion of 
ortho- and para-H 2 as well as non-LTE line radiative transfer, 
we model the observed DCO^ emission. We find that using a 
constant CO abundance predicts a DCO^ peak at around 4", 
twice further out than observed, irrespective of the adopted o/p 
ratio and density profile. Our model results show that a drop CO 
abundance with a decreased temperature profile by a factor of 
1.5 generates a peak at the same position as the observed emis¬ 
sion. Thus, the observed DCO'^ peak is closer to VLA1623A 
and with a lower temperature (11-16 K) than that expected from 
a spherically symmetric physical structure constrained by con¬ 
tinuum data and source SED. 

3. The observed DCO^ and C'^O emission are not described 
by the same physical structure. In our model, the C**0 emission 
is well reproduced by a drop CO abundance with the density 
profile increased by one order of magnitude at 1" (<120 AU) 
radii, in comparison with the density profile needed to reproduce 
the observed DCO^ emission. A constant CO abundance and a 
decreased temperature profile over-predict the extent of the C'^O 
emission. 

4. Disk-shielding is the best possible explanation for the 
observed DCO^ emission toward VLA1623A. Disk-shielding 
causes the inward shift of the CO freeze-out region along the 
plane of the disk, lowering the dust temperature to <20 K, gen¬ 


erating a ring of molecules whose abundance is enhanced by 
low temperatures such as DCO^. The rest of the envelope is 
largely unaffected by the disk, thus the CO freeze-out shell pre¬ 
dicted by spherically symmetric radiative models is expected to 
be located further out. This prediction is confirmed by our recent 
ALMA Cycle 2 ACA observations, which show that the DCO'^ 
emission along the outflow axis lies at larger radii, ~5", consis¬ 
tent with constant CO abundance models without any alteration 
to J0rgensen et al. (2002|l’s temperature and density profile of 


VLA1623. 

5. The disk-envelope interface in VLA1623A is shown to 
have a broken transition in density and temperature, with the im¬ 
pact generated by the presence of a disk being observable from 
small to large scales. 

Our observations and modeling results for VLA1623A show 
the disk-envelope interface to have different physical conditions 
than other parts of the envelope. Our results also highlight the 
drastic impact that the disk has on the temperature structure 
at ~100 AU along the plane of the disk, with the effect be¬ 
ing observable even at large scales. We suggest performing fur¬ 
ther observations to determine whether the unequal physical and 
chemical conditions observed in the disk-envelope interface of 
VLA1623A is a common phenomenon in protostellar systems 
with rotationally supported disks or a special condition of the 
present source. 
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Appendix A: Modeling results: additional figures 

Here we present additional figures illustrating details discussed 
in Section [42| and |4.3| Figures [AT| and [A!^ show the effect of 
CO abundance Xco in the drop on the DCO^ peak for Test 1 and 
4, respectively. Figure [A3] highlights the effect of the upper and 
lower o/p limits on the DCO^ abundance for the best approxima¬ 
tion models to the DCO^ observations. The effect of the factor 



different tests, CO abundance prohles and cases. 
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Fig. A.l. Chemical network results for Drop CO abundance with the Test 1 source profile from |j0rgensen et al.| ( |2002| l for four cases (Table [^. 
Case d is the same as in the top row center panel of Fig.|^ 



Fig. A.2. Chemical network results for Drop CO abundance with the decreased-temperature (Test 4) source profile for four cases (Table|^. Case 
d is the same as in the bottom row center panel of Fig.|^ 


Drop: Case d, upper o/p Drop: Case d, thermal o/p Drop: Case d, lower o/p 




Fig. A.3. Chemical network results for Drop CO abundance with the decreased-temperature source profile (Test 4) for case a with upper-limit 
(left), thermal {center) and lower-limit {right) o/p ratio. 
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Fig. A.4. DCO^ Cycle 0 observations compared with the radial profile of the constant and drop CO abundance models, showing the effect of the 
factor used to decrease the temperature, for case d. A factor of 1.5 with drop CO abundance provides the best approximation to the observed radial 
profile. 
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Fig. A.5. Comparison of the observations and radiative transfer modeling of the chemical network results. The black line shows the DCO^ 
observations integrated over the southern clump, with the gray solid lines showing the Icr and 4cr levels, respectively. Models are referred to by 
test, abundance profile (Constant "C" or Drop "D") and case. Color lines show tests 1 (blue), 2 (green), 3 (orange) and 4 (purple) for the top left 
panel, and cases a (blue), b (green), d(orange) and f (purple) for the other four panels. 
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Fig. A.6. As in Fig. 


A.5 


but for C'*0. Gray solid and dashed lines show the Icr and 3cr levels respectively. Note that constant CO abundance 
over-predicts the amount of C**0 emission. 









































